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The results of a high-statistics study of scalar and pseudoscalar meson prop- 
agators in quenched lattice QCD are presented. For two values of lattice spac- 
^ • ing, P = 5.7 (a ~ .18 fm) and 5.9 (a ~ .12 fm), we probe the light quark 



mass region using clover improved Wilson fermions with the MQA pole-shifting 
ansatz to treat the exceptional configuration problem. The quenched chiral 
loop parameters tjiq and are determined from a study of the pseudoscalar 
hairpin correlator. Prom a global fit to the meson correlators, estimates are ob- 
tained for the relevant chiral Lagrangian parameters, including the Leutwyler 
parameters L5 and Lg- Using the parameters obtained from the singlet and 
nonsinglet pseudoscalar correlators, the quenched chiral loop effect in the non- 
singlet scalar meson correlator is studied. By removing this QCL effect from 
the lattice correlator, we obtain the mass and decay constant of the ground 
state scalar, isovector meson oq. 
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1 Introduction 



Improved methods for studying the regime of smah quark mass in lattice QCD provide 
the reaUstic prospect of quantitatively determining the parameters of the low energy chiral 
Lagrangian of QCD from first principles. Although a definitive comparison with experi- 
ment requires the analysis of full and/or partially quenched simulations, detailed studies 
of chiral behavior in the quenched approximation are of interest for several reasons. First, 
the characteristic quenched chiral loop effects which arise from the anomalous double-pole 
structure of the quenched, flavor-singlet pseudoscalar propagator will also occur in par- 
tially quenched calculations (at a level determined by the mismatch between valence and 
sea quark masses) pp. The observation of these anomalous effects in the quenched theory 
should provide a useful baseline for future chiral analysis of full QCD. Second, although it 
is not a unitary theory, the quenched approximation can be analyzed in an effective La- 
grangian framework 21 El) yielding a well-defined set of low-energy constants in quenched 
chiral perturbation theory. (In practice, this requires the assumption that the Ua{^) break- 
ing from the anomaly can also be treated perturbatively). Comparison of these constants 
with those of full QCD can provide valuable insight into the role of closed quark loops 
in hadron phenomenology. In addition, the study of chiral behavior in the quenched ap- 
proximation provides useful information about the interplay between topological charge 
and chiral symmetry breaking in QCD. For example, the Witten-Veneziano formula relates 
the gluonic component of the rj' mass in full QCD to the topological susceptibility of the 
quenched theory. 

In two previous papers [IJ [5] we reported results of a study of the chiral behavior 
of scalar and pseudoscalar meson propagators in quenched QCD at /? = 5.7 using clover 
improved Wilson fermions. In this paper we present results from a new data set at /? = 
5.9 in], compare them with the /3 = 5.7 results, and summarize the main conclusions of 
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this study. We also compare our results to those of other recent studies. [7j An important 
ingredient in our analysis is the use of the MQA pole-shifting procedure [HI to resolve 
the exceptional configuration problem. Our experience with this technique has led us to 
conclude that it provides a practical and quantitatively acceptable resolution of the problem, 
which eliminates the spurious statistical fluctuations of exceptional configurations without 
systematically biasing the final results. This conclusion is based on the consistency between 
the light-quark results and those of heavier quarks where the pole-shifting has a negligible 
effect on the propagators. It is also supported by the overall agreement we observe between 
our results (in and out of the pole region) and theoretical expectations based on quenched 
chiral perturbation theory. (See Section 6) 

Since the introduction of the MQA procedure, other methods for avoiding the excep- 
tional configuration problem have been explored. These include twisted-mass QCD jlUj and 
the use of exactly chiral (overlap or domain- wall fermions. All these approaches 
have in common the fact that Wilson-Dirac eigenvalues at positive real quark mass are 
eliminated, thus resolving the problem. It should be noted that exactly real eigenmodes of 
the Wilson-Dirac operator, which are the cause of the exceptional configuration problem, 
make a negligible contribution to physical quantities in the infinite volume limit (vanishing 
like Thus any prescription which effectively removes these poles from the physical 

region should provide a satisfactory resolution of the problem for sufficiently large volume. 
The MQA pole shifting is a minimal prescription for accomplishing this. A more stringent 
test of the procedure is the study of chiral behavior for very light quarks in a finite volume 
which is large compared to the QCD scale but comparable to the chiral scale. In this regime, 
finite volume effects are large but calculable in QxPT, simply by replacing loop integrals 
by finite- volume momentum sums. As discussed in the scalar, isovector (valence) meson 
propagator exhibits a prominent quenched chiral loop effect arising from the tj'-tt interme- 
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diate state. For the lightest quark masses we study, the finite volume effects expected from 
QxPT are quite large. Thus, the detailed agreement (as a function of both time and pion 
mass) between the measured scalar propagator and the finite-volume one-loop calculation 
provides a convincing demonstration that the MQA procedure is an effective method for 
exploring the light quark regime with Wilson fermions. 



2 Quenched Chiral Perturbation Theory 

To analyze the quenched theory in a chiral Lagrangian framework, one introduces wrong- 
statistics ghost quark fields to cancel closed loops, yielding a low-energy chiral Lagrangian 
with a graded C/(3|3) x [/(3|3) symmetry (5]. At the one-loop level, this is equivalent to 
the simpler and more direct approach to quenched j^j which begins with an ordinary 

U{3) X U (3) chiral Lagrangian describing a nonet of Goldstone bosons. To leading order, 
this is 



Trid^U^df^U) + Tr{x^U + Wx)] (1) 



where ?7 is a U{3) x U{3) chiral field and x is the pseudoscalar mass matrix. Our analysis 
also incorporates the following fourth-order terms in the chiral Lagrangian jl4j . 

£4 = L^Trid^U^d^Uix^U + Wx)) + LsTr{x^Ux^U + WxUh) (2) 

The effect of the axial U(l) anomaly is introduced as an explicit symmetry breaking term 
consisting of a flavor-singlet pseudoscalar rj' mass term and a field renormalization, 

^hp = ^(«<l.5V5^7?'-m2r?'2) (3) 



where 

^ = 2 



Tr ln([/t) - iTr ln([/)) (4) 
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Finally, we will also analyze the scalar, isovector meson propagator, which turns out to be 
well-described by a combination of a heavy oq meson and an r/'-vr loop diagram. Thus, we in- 
corporate a scalar-isovector meson field, using the formalism of nonlinear chiral Lagrangians 

C,c = ^tr{DaDa] - Knltr{aa] + fstr{x^VUaVU + xVU^(^VU^} (5) 

where D is a chirally covariant derivative. One of our motivations for studying the scalar 
correlator is the expectation of a prominent quenched chiral loop effect from the ij'-ir inter- 
mediate state, as discussed in Ref. jSj. The agreement between the lattice correlator and 
the one-loop calculation is very good, particularly for the (3 = 5.9 results, as discussed in 
Section 6. 

To summarize, the low energy chiral Lagrangian used in our analysis is 

£ = C2 + C.A + Chp + Csc (6) 

In the quenched approximation, we have the supplementary rule that multiple r/' mass 
insertions on a given pseudoscalar line are excluded. More generally, any xP'^ diagram 
corresponding to a quark-line diagram with internal closed loops is discarded [3]. At the 
one-chiral-loop level, these rules are unambiguous, and equivalent to the more systematic 
procedure of introducing ghost fields P]. 

3 Lattice Parameters 

The calculations discussed in this paper were carried out on the Fermilab ACPMAPS and on 
the UVA Linux cluster GARCIA. The two Monte Carlo gauge ensembles analyzed consisted 
of 300 configurations at /? = 5.7 on a 12^ x 24 lattice, and 350 configurations at /? = 5.9 on a 
16'^ X 32 lattice (the Fermilab b and c ensembles). Quark propagators were calculated with 
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clover improved Wilson action. The clover coefficients used were Cgw = 1-57 for /3 = 5.7 
and Csw = 1-50 for (3 = 5.9. For (3 = 5.7, the quark propagators were calculated for 
K = .1410,.1415, .1420, .1423, .1425, .1427, and .1428, with = .14329, while for /? = 5.9, 
propagators were calculated with k = .1382, .1385, .1388, .1391, .1394, and .1397, with Kc = 
.14013. In physical units, this corresponds to a range of pion masses of 275 to 565 MeV 
for P = 5.7 and 330 to 665 MeV for (3 = 5.9. Here and elsewhere, we will quote results in 
physical units using the rho mass to set the scale. 

An analysis of smeared and local rho propagators on our ensembles yields mpa = .690(8) 
and .469(3) for /3 = 5.7 and f3 = 5.9, respectively. An analysis of the axial-vector meson 
channel also yields a mass for the oi meson of 1.15(7) and 0.77(3) for (3 = 5.7 and (3 = 5.9 
respectively. Using the rho mass to fix the scale gives = 1.12 GeV for (3 = 5.7 and 1.64 
GeV for (3 = 5.9. The resulting physical mass for the oi (1290 MeV for f3 = 5.7 and 1260 
MeV for (3 = 5.9) is close to the mass of the observed ai(1260) resonance. 

To get some idea of the systematic error associated with choice of scale, we will some- 
times quote equivalent results using the charmonium IS- IP splitting scales of 1.18 GeV and 
1.80 GeV for the two ensembles. The MQA pole-shifting procedure [S] was applied to all 
quark propagators. For (3 = 5.7, all poles below k = .1431 were located and shifted, while 
for (3 = 5.9, all poles below k, = .1400 were shifted. 

4 The Hairpin Insertion, rj' Mass, and Topological 
Susceptibility 

We begin by determining the parameters uiq and a<j> in the term Chp, Eq. (j^J. These 
parameters are extracted from the two-quark-loop ("disconnected") piece of the quenched 
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flavor singlet pseudoscalar correlator, 

Ah{x) = {Tn'^Gix, x)Trj^G{0, 0)) (7) 

The lowest order QxPT approximation to ^ is the tree graph with a single hairpin insertion 
between two pion propagators. In momentum space, this is 



= fP ^2^.^2 ("^0 + . ^ o fP (8) 



1 / 2 2\ 1 

pZ _|_ V / _|_ j^z 

where fp is the pseudoscalar decay constant, 

fp = (OlV^T'V'k) (9) 



Fourier transforming over pQ and setting ^ = 0, we have 

Ami 



Ah{p = 0,t) = -f^ [C+ + C-m^t] e-"*-* + {t ^ T - t) (10) 



where 

C± = ml±a,^ml (11) 

In our previous analysis of the /? = 5.7 ensemble, the hairpin correlator was studied 
for both local and smeared sources and compared with the pion pole residues of the corre- 
sponding valence propagators. This analysis demonstrated a remarkable absence of excited 
state contamination in the hairpin correlator, even when the sources were only separated 
by one or two time slices. Moreover, the time-dependence of Ah{t) was well-described at 
all times t > 2 by the formula ()1U() with a$ = 0, i.e. by a pure momentum-independent 
mass insertion. For the present analysis, we have fit both the 5.7 and 5.9 ensembles to the 
full two-parameter formula (|1U() in order to obtain an accurate estimate of a^. For the 5.7 
ensemble, acceptable x^'s were obtained by fitting a range of times from t = 3 to 12. Using 
the fully correlated error matrix, the covariant ^or these fits ranged from 0.8 to 1.5 per 
degree of freedom. For the hairpin correlators at /3 = 5.9, we obtained very good fits to the 
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formula H10|) over the entire time range from t = 1 to 16. Here the correlated x^'s ranged 
from 0.3 to 0.5 per degree of freedom. An example of a hairpin fit for [3 = 5.9 and n = .1394 
is shown in Fig. ^ The solid line is the pure dipole fit with q$ = 0. In Tables Q ^ind 
l^we give the /3 = 5.7 and 5.9 results for mo and 0$. Also shown in the last column are 
the values of obtained from the 1-parameter pure dipole fit with 0$ = 0. Considering 
first the results for a^, the values for the 5.7 ensemble are negative by about one to two 
standard deviations, while the values for (3 = 5.9 are slightly positive, also by about two 
standard deviations. The values for different k's within each ensemble are highly correlated, 
so the deviation of a<j> from zero in either data set has little statistical significance. Ignoring 
K-dependence and averaging the values within each ensemble, we get 

acjb = -0.15 ± 0.10, /3 = 5.7 (12) 

and 

«$ = 0.05 ± 0.03, /3 = 5.9. (13) 

If we ignore any possible lattice spacing dependence and average the two data sets in 
quadrature, we get the final result 

a$ = 0.03 ± 0.03 (14) 

This can be regarded as a success of the large- A'^c view of the anomaly where this renormal- 
ization is an order l/N^ effect. In the subsequent analysis we will take 0$ = 0. 

Fitting the hairpin correlators to the pure dipole = form (last column of Tables ^ 
and 121) we obtain the chirally extrapolated values (in lattice units) 

mo = .348(4), /? = 5.7 (15) 

and 

mo = .232(4), (3 = 5.9. (16) 
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Table 1: Fit parameters mo and a<i> for the /3 = 5.7 hairpin correlators. All masses 
are in lattice units. 





m^ 


mo 


a$ 


mo(a;$ = 0) 


.1410 


.505(2) 


.269(26) 


-.17(10) 


.280(10) 


.1415 


.450(3) 


.291(24) 


-.18(10) 


.294(10) 


.1420 


.386(3) 


.310(21) 


-.19(10) 


.308(10) 


.1423 


.342(4) 


.321(19) 


-.19(10) 


.316(10) 


.1425 


.307(4) 


.326(19) 


-.16(11) 


.321(10) 


.1427 


.267(5) 


.326(19) 


-.10(12) 


.322(11) 


.1428 


.245(6) 


.323(19) 


-.03(13) 


.322(11) 



Table 2: Fit parameters mo and for the (3 = 5.9 hairpin correlators. All masses 
are in lattice units. 





m^ 


mo 


Q;$ 


mo(a$ = 0) 


.1382 


.411(3) 


.188(14) 


.04(2) 


.194(5) 


.1385 


.378(3) 


.192(13) 


.04(2) 


.198(5) 


.1388 


.343(4) 


.197(11) 


.04(2) 


.203(5) 


.1391 


.304(4) 


.204(11) 


.05(2) 


.209(5) 


.1394 


.261(5) 


.211(9) 


.06(3) 


.217(5) 


.1397 


.204(5) 


.220(9) 


.11(4) 


.226(6) 
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Using the rho mass scale and including a flavor factor of \/3, this gives the gluonic component 
of the rl mass 

= 659(12) MeV, /3 = 5.9 (18) 
If we instead use the charmonium scale we get 

^giue ^ 7^2(9) MeV, f] = 5.7 (19) 
= 723(13) MeV, (3 = 5.9 (20) 

We conclude that the rj' mass scales reasonably well between /3 = 5.7 and 5.9, well within 
the systematic uncertainty associated with different ways of determining the lattice spac- 
ing. The values we obtain for the rj' mass insertion are somewhat low compared to the 
estimate of ~ 850 MeV obtained from the physical ij' mass and chiral perturbation the- 
ory. Although we see approximate scaling, it would require calculations at larger values 
of P to rule out a significant lattice spacing effect. It is also worth remembering that the 
whole framework in which the quenched hairpin diagram is interpreted as a mass insertion 
is only demonstrably valid in the limit of large Nc, so some discrepency between the lattice 
calculation of rriQ and the phenomenological estimate might be expected. A recent calcu- 
lation of the r]' mass in two-flavor full QCD by the CPPACS collaboration ^3] gave the 
result m^/ = 960(87)1^286^^^' good agreement with experiment. Detailed comparisons 
between quenched and full QCD studies of the r]' should provide a better understanding of 
the accuracy of laige-Nc arguments in the framework of chiral Lagrangians. 

The overall size of quenched chiral loop effects is determined by the parameter 5, which 
can be computed from the hairpin insertion mass mo and the axial vector decay constant 
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/a (evaluated in the next Section), 



5 




(21) 



247r2/2 



Using the chirally extrapolated values of mo and /a, we obtain 



5 = .099(3), P = 5.7 



(22) 



and 



5 = .108(4), ^ = 5.9. 



(23) 



It is interesting to consider not only the value of 5 in the chiral limit, but also the effective 
value of 5 at a given quark mass by computing the quantity ()21() from the values of tjiq 



plotted in Fig. 2. This plot shows a rather strong quark mass dependence of the effective 
QCL parameter, which may provide at least a partial explanation of the fact that many 
of the determinations of 6 from lattice studies of quenched chiral logs [7] have favored a 
value of 5 substantially smaller than the phenomenological estimate of 5 ~ 0.17. From Fig. 
1^1 we see that, for pion masses > 300 MeV where most studies have been carried out, the 
value of 6eff is smaller than the value in the chiral limit by as much as a factor two. The 
decrease of Jg// with increasing quark mass represents the combined effect of a decreasing 
value of niQ and an increasing value of /a as the quark mass increases. Although the 
negative slope of 6eff has the effect of suppressing quenched chiral logs, it is nevertheless 
more consistent to treat 5 as a constant in fitting to chiral Lagrangian parameters, since 
the effective mass dependence should arise from higher order terms in the chiral expansion. 
This is the procedure we adopt in the subsequent analysis of the scalar and pseudoscalar 
correlators, where the best QxPT fit favors a value of 5 about half as large as that obtained 
from the chirally extrapolated hairpin result. We might expect to find a larger value of 5 



and /a at that mass. The values of (5e// vs. pion mass' 



? for both (3 = 5.7 and 5.9 are 
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if studies were carried out well below m-,^ = 300 MeV. [It is interesting that a recent study 
jl5j using overlap fermions, which went as low as m-,^ = 180 MeV, found a large value of 6. 
However, the value 6 = 0.26(3) obtained in Ref. much larger than even our chirally 

extrapolated result of 0.108(4), indicating that there are other systematic differences in the 
calculations. Further chiral studies comparing different fermion actions on the same gauge 
configurations would be of considerable interest.] 

The "allsource" quark propagators used to calculate the hairpin correlators jl6j can also 
be used to calculate the topological susceptibility. Using the integrated anomaly method 
[111 13 , we calculate a winding number for each gauge configuration from the pseudoscalar 
charge integrated over the whole lattice. From these winding numbers, we compute the 
topological susceptibility xt = (j^^)/^- Using the rho scale, this gives 

Xt = (178(4) MeV)^ /3 = 5.7 (24) 
= (171(3) MeV)^ p = 5.9. (25) 

Results quoted previously ^ used the charmonium scale, which gives 

Xt = (188(4) MeV)^ /3 = 5.7 (26) 
= (190(3) MeV)^ p = 5.9. (27) 

5 Pseudoscalar Masses and Decay Constants 

The one-loop chiral Lagrangian analysis of the pseudoscalar and axial-vector propagators 
for the (3 = 5.7 ensemble has been described previously @]. Here we briefly review that 
analysis and compare the previously reported results with the new results at /? = 5.9. 
We also compare with results of the Alpha collaboration |S] and discuss some issues 
associated with extracting the Leutwyler parameters L5 and Lg, or correspondingly 05 and 
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Figure 2: The quenched chiral log parameter Sejf vs. squared pion mass. Results 
from both (3 = 5.7 (x's) and 5.9 (boxes) are plotted. The linear fits include all mass 
values for j3 = 5.7 and the four lightest masses for /5 = 5.9. 
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as- (Note: Here and elsewhere, we use a conventional notation for the rescaled parameters: 
ai = 8(47r)2Li). 

For both ensembles, we calculated propagators using smeared pseudoscalar, local pseu- 
doscalar, and local axial-vector sources and sinks. The calculations were done for all meson 
propagators with both degenerate and nondegenerate quark masses. The chiral Lagrangian 
parameters were extracted from a global fit to all pseudoscalar masses and decay constants 
based on one-loop quenched xPT for the Lagrangian £2 + ^^4-1- Chp, as discussed in Section 
3. For the lightest pion masses we studied, finite volume effects on chiral loop integrals are 
potentially significant, so all one- loop calculations were carried out with the appropriate 
finite-volume momentum sums rather than loop integrals. In addition, quadratic and log- 
arithmically divergent integrals are regularized by subtraction at a cutoff scale A ~ To 
summarize, a generic loop integral of the form 

^ r,, I 1 



is replaced by a cutoff momentum sum 

= 167r2 ^ (d{p,, Mi)D{pj,M,) - D{p, Af) (29) 
p 

while a quadratically divergent integral 

Ii = \ I (30) 

is replaced by 

li = levr^ ^ (^D{p, Mi) - D{p, A) - (A^ - M^)D{p, Af) (31) 
p 

In these expressions, D{p,M) is the free boson propagator and the momentum sums are 
defined according to the physical size of the corresponding lattice volume. With these 



15 



0.20 



0.15 



0.10 



0.05 



0.00 



J I I I I I I I I I I I 1_ 



I . . 



_l I 1_ 



0.00 0.01 0.02 0.03 0.04 0.05 0.06 

quark mass 



Figure 3: Pion mass squared for equal quark masses for the the (5 = 5.9 ensemble. 

definitions, the value of the squared pseudoscalar meson mass up to first order in L5, Lg, 
and 5 is 

1 



+ 6{iij + lij)] + j^SL^XijSJij} 



(32) 



where tq is a slope parameter, Xi = 2romj, Xij = {Xi + Xj)/2 and m,/ = ln(l + l/(2K;j) 
1/(2k;c))- Here lij = {luXi + IjjXj)/Xij, 



J,, = (Ji + Ij - {Mt, + M^j)hj) /2, (33) 

Jij = Jij/Xij- The loop integrals lij are defined in 0]. The fits for the pseudoscalar masses 
are shown in Figures |30 

For the pseudoscalar decay constants: 

/P;,,= y2/ro(l + 0.255(/ii + Ljj + 2Ii,)){l 
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Figure 4: Chiral slope parameter scan for /5 = 5.9. Individual chiral slopes are 
denoted by (x). The scan number is 6 * (i — 1) + j for the pair of k, values (kj, kj). 
The six k, values are ordered from lightest to heaviest quark mass. Also shown are 
the ratios to the chiral log fit (boxes) which yields a global average of 1.712(39). 
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+ j^{4Ls - L5)X^A1 + SiiiJ + lij)] 

- ^L^Sx^A{k^ + hj) - {Ju + 4)]} (34) 

The fits for the pseudoscalar decay constants are shown in Figures I5I6I 
For the axial vector decay constants, we have: 



4 

7 



{ 1 + -p^XiM^ + ^{'^hj - iij)]} (35) 



The fits for the axialvector decay constants are shown in Figures I7I8I In Figures ^ El and 
IHlwe show both the measured values and the ratio values where the fitted chiral log factors 
in Eqs. WVA IIMI and have been divided out. The results are scanned over all combinations 
of the quark mass values. 

The chiral Lagrangian parameters listed in Table 01 are the result of a global correlated 
fit of the xPT^ expressions to masses and decay constants with both equal and unequal quark 
masses. The results obtained from this global fit are generally consistent with evaluations 
extracted from more limited fits to equal quark mass data. In particular, the value of L5 
may be estimated directly from the axial vector decay constant For equal quark masses 
/a has no quenched chiral logs, so the mass dependence should be linear in the chiral limit, 

fA = A + Bml (36) 

(Note that, in our notation, includes a factor \/2 relative to /jr? i-e. it's physical value is 
132 MeV.) A determination of -L5 is obtained from the product of slope and intercept: 

^5 = ^ (37) 

Note that A has units of mass and B has units of {mass)~^, so that L5 is dimensionless 
and can be evaluated without reference to a mass scale. We obtain the results (in GeV 
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Figure 5: Pseudoscalar decay constant, fps, versus pion mass squared, for the 
(3 = 5.9 ensemble. Only values for equal quark masses are shown. 
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Figure 6: Pseudoscalar decay constant scan for f] = 5.9. Also shown are the ratios to 
the chiral log fit which yields a global average fps = 0.1501(77). The values 
are ordered as in FiglH 
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Figure 7: Axial vector decay constant, /ax, versus pion mass squared, for the 
/3 = 5.9 ensemble. Only values for equal quark masses are shown. 
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Figure 8: Axial vector decay constant scan, /ax for = 5.9. Also shown are the ratios 
to the chiral log fit which yields a global average /ax = 0.0915(13). The (kj, Kj) values 
are ordered as in FiglH 
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using the rho scale) 

/a = (.1653(17) + .m{5)ml) x Za GeV, /3 = 5.7, (38) 

and 

fA = (.1515(20) + .106(7)m^) x Za GeV, f3 = 5.9. (39) 

These fits are shown in Figs. |ll|Hland|Sl Using Za = 0.845 and 0.865 for /3 = 5.7 and 5.9, 
respectively, we find 

L5(103) = 1.72(11), 05 = 2.18(14), /3 = 5.7 (40) 

and 

L^iW^) = 1.50(9), as = 1.89(11), f3 = 5.9 (41) 

These values are consistent with the global fit results in Table |5J 

A recent analysis of quenched data from the Alpha collaboration|Sj obtained a value of 
cts = 0.99(6). Their analysis of the /a ratios, Rf{x), neglects a factor of 

, (42) 

(1 +yref2«5) 

(c.f Eq. (3.7)-(3.9) of [S]) reflecting the difference between xPT expansion about zero mass 
and their use of a large reference mass. Including this factor increases their estimate of as 
to 

as = 1.18(7). (43) 

If only statistical errors are considered, this estimate is still significantly below our values. 

We can identify a number of differences in the two analyses which could account for 
much of the discrepancy. First, our data uses somewhat coarser lattices and an independent 
determination of the clover coefficient. This could affect the comparison with the Alpha 
analysis although they see little lattice spacing dependence in their data. There are also 
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Table 3: Best fit for the chiral Lagrangian parameters for the (3 = 5.7 and f3 = 5.9 
lattices. 



Parameter 


p = 5.9 


(3 = 5.7 


f=fn 


0.091(2) 


0.100(2) 


Ls X 10^ 


1.55(27) 


1.78(35) 


(Lg - L5/2) X 103 


0.02(5) 


0.14(7) 


ro 


1.71(8) 


1.99(12) 


6 


0.053(13) 


0.059(15) 



differences in methodology. Our results are obtained from a linear fit to /_4(m^) where the 
value of L5 is determined from the dimensionless product of the slope and the intercept 
of the fit. This method is insensitive to physical scale parameter but is quite sensitive to 
the value of the axial-vector renormalization constant, being proportional to Z^. In the 
context of the tadpole renormalization scheme ^Sj our results in Eqs. I4UI41I include the 
perturbative renormalization constants Za_ = 0.845 and Za_ = 0.865 at /? = 5.7 and (3 = 5.9, 
respectively. These values are derived from the perturbative formulas given in Ref. |19j . 

As an alternative method, we could determine the renormalization factors by requiring 
a fit to the physical value of /yr = 93 MeV. The result will now be sensitive to the choice 
of scale parameter. Fixing the scale with the rho mass, the renormalization constants are 
given by Za = 0.80 for (3 = 5.7 and Z^ = 0.87 for (3 = 5.9. Using these renormalization 
constants modifies our predictions to 

L5(103) = 1.54(10) , 05 = 1.95(13), (3 = 5.7 (44) 
^5(10^) = 1.52(9) , 05 = 1.91(11), (3 = 5.9 (45) 

This result is still substantially above the value obtained in the Alpha analysis. We note 
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that, in our results, the mass dependence of the tadpole renormalization factor contributes 
about 17% to the slope of /^(m^). For linear extrapolations of /yi(m^) the above method 
is actually exactly equivalent to the Alpha ratio method with the inclusion of the correction 
factor (Eq. H^ . 

In the Alpha procedure, the physical value of /tt is used along with a particular choice 
of scale parameter. (This procedure, which involves the physical value of /^r, and/or nip, is 
somewhat adhoc since the quenched values for f^^ and nip need not agree precisely with the 
Particle Data Group.) We note that caution is required in any direct comparison between 
the quenched and unquenched values of as . We have emphasized that there are no quenched 
chiral logs which affect the extrapolation of /^(jn^) for equal quark masses. However there 
is a strong chiral log effect in the unquenched theory. The Gasser and Leutwyler analvsisjl4j 
of the one-loop chiral logs determines 05 from the physical values of / fn 

as = 2.8 ±0.6-3 ln(^f/"^»?)- (46) 

Hence, as = 0.5 for a large scale, fi = 471/7^, 05 = .99 for /i = 1 GeV and as = 1.76 
for = nip. The quenched theory really corresponds to the leading order of the l/Nc 
expansion of the full theory which can not be isolated phenomenologically from the higher 
order contributions. 

6 The Quenched Scalar Propagator 

One of the most dramatic effects of quenched chiral loops is observed in the scalar valence 
propagator. The behavior of this propagator for the [3 = 5.7 ensemble was studied in jH]. 
There it was found that the propagator was well-described as a combination of a short-range 
positive exponential associated with a heavy (> 1 GeV) scalar qq meson state and a long- 
range negative tail arising from the ?/-7r intermediate state. In the quenched approximation, 
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the r/'-TT loop diagram exhibits not just a quenched chiral logarithm, but a quenched chiral 
power, with infrared behavior ~ d'^p/p^. This long-range component is well-described in 
both shape and magnitude by the finite-volume i]'-tt loop calculation. [Note: In the /? = 5.7 
results, the ~ 2 standard deviation discrepency between the data and the calculation 
for t > 7 (c.f. Fig. 10 of Ref. jHI) we now believe to be a statistical fluctuation. See 
below.] As discussed in jS] the one-loop term that dominates the scalar propagator at 
small quark mass is determined with no adjustable parameters by the chiral Lagrangian 
parameters 771,^,/, and tq, already fixed from the analysis of the pseudoscalar sector. In 
particular, the long-range scalar propagator exhibits a very strong mass dependence in the 
light quark regime, which is very well explained by the dependence of the finite volume one 
loop contribution on in^. 

The analysis of the scalar correlator for /3 = 5.9 confirms the main conclusions of Ref. 
(See Figures f9ll0() For the 5.9 ensemble, we find excellent agreement between the long-range 
behavior of the correlator and the QxPT loop calculation (which contains no adjustable 
parameters), with the theoretical formula matching the data all the way out to the largest 
time separation of t = 16. The details of the analysis of the (3 = 5.7 data were discussed in 
[3]. The f3 = 5.9 ensemble was analyzed similarly, though with some simplifications. First, 
in [Sj both smeared and local sources were used, and a factorization procedure allowed 
an estimate of the mass of the excited Oq meson. This produced a good fit to the scalar 
correlator all the way down to t=l. For the analysis of the 5.9 data, we have not carried 
out an independent estimate of the excited Oq contribution, but have simply rescaled the 
parameters obtained from the 5.7 analysis, holding this contribution fixed in the fits which 
determine the ground state cq parameters. For the 5.9 data we used a time range t > 3, 
so the fits were insensitive to the choice of excited state parameters. The other difference 
in the 5.9 analysis is that the lightest pion mass (330 MeV) is somewhat heavier than that 
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Figure 9: Isovector scalar correlator for the (3 = 5.9 16^ x 32 lattices. The lightest 
.1397) (red o) and heaviest ( .1382) (blue x) correlators are 

shown. 
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for the 5.7 study (275 MeV). In the formula used to parametrize the rj-ix contribution 
was obtained by resumming iterated bubble graphs to all orders. In the analysis of the 5.9 
data, we have found that only the one-loop graph is significant, so we have discarded higher 
order bubbles in the formula used to fit the correlator. Thus, defining the time-dependent 
zero-momentum scalar correlator as 

A(t)^X](V5iV2(^,t)Vi2V'i(0)) (47) 

X 

we fit the lattice correlator to the function (c.f. Eq. (14) of Ref. [Hj) 



A(i) ~ 32r2-^e-'-^* + ^rlBnv{t) (48) 



where -B(t) is the p = Fourier transform of the rj-ix bubble graph, calculated in a finite 
volume: 

. , 1 ^ — ^ 1 TTIj^ . , 

Bh^ [p) = —^ ______ + )2 (49) 

Note that 771,^, tq, and mo are already determined from the pseudoscalar correlator analysis, 
so the only two fit parameters are the scalar mass and decay constant nis and fs- In Fig. 
we show the data for the scalar correlator at (3 = 5.9 and the lightest quark mass, k, = .1397, 
along with a fit given by the function (|48|) . 

The mass of the ao meson obtained from our fits is plotted in Fig. 1111 In the chiral 
limit we obtain (using the rho scale) 

niao = 1285(60) MeV at /? = 5.7 (50) 
= 1326(86) MeV at (3 = 5.9 (51) 



In the chiral loop calculation of the r/'-vr intermediate state contribution, we have in- 
cluded the effect of finite physical volume, replacing momentum integrals by discrete sums. 
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Figure 10: The scalar correlator for /3 = 5.9 and k = .1397. The solid line is a fit 
consisting of the sum of a scalar meson pole and an ri'-n loop diagram. The dashed 
line is the loop term in the infinite volume limit. 
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Figure 11: Mass of the Oq meson for f3 = 5.7 (boxes) and f3 = 5.9 (x's). 
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Figure 12: The scalar decay constant of the Oq meson for /? = 5.7 (boxes) and /? = 5.9 
(x's). 

For the hghtest pion masses the finite volume effect is quite large, due to the infrared be- 
havior of the loop contribution. It is interesting to note that the Monte Carlo results for the 
scalar correlator agree quite nicely with the prediction of finite volume QxPT, but disagree 
significantly with the corresponding infinite volume calculation. This strong finite volume 
effect is also seen by the RBC collaborationpU] . The dashed line in Fig. ^Jis the infinite 
volume loop calculation. For t > 5, the correlator is dominated by the r^'-vr loop, and it 
is clear that the Monte Carlo result exhibits the predicted enhancement from the finite 
volume effect. From the point of view of Dirac eigenmodes, the finite volume dependence 
predicted by xPT arises from a subtle interplay between exact zero modes and near zero 
modes. In the infinite volume limit, the contribution of exact zero modes vanishes, but for 
finite volume these modes are essential for reproducing the correct chiral behavior. Since 
the MQA pole shifting ansatz consists of repositioning some exact zero modes, it is reassur- 
ing to see that the resulting correlators exhibit excellent agreement with the finite volume 
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effect predicted by QxPT. 

Finally, the scalar decay constant, fs, of the meson obtained from our fits is plotted 
in Fig. 1121 In the chiral limit we obtain (using the rho scale) 

= 64(5) MeV at (3 = 5.7 (52) 
= 68(3) MeV at (3 = 5.9 (53) 

7 Conclusions 

The lattice calculations described in this paper and our previous works [^[3115] have focused 
on the chiral properties of meson correlators in the flavor singlet and nonsinglet pseudoscalar 
and nonsinglet scalar channels. In addition to exhibiting for the first time a number of 
anomalous chiral effects due to quenching, the results have confirmed a level of overall 
consistency of the low energy chiral Lagrangian description of meson properties in QCD, 
and provided quantitative estimates of the relevant chiral Lagrangian parameters. 

A central feature of this study is the accurate calculation of the pseudoscalar flavor- 
singlet hairpin insertion responsible for the gluonic component of the r/' mass. In this 
calculation the MQA pole shifting ansatz has a particularly salutary effect. Since the 
pseudoscalar hairpin insertion arises from the Ua{^) anomaly, it is particularly sensitive to 
topological features of the gauge field. As a result, the exceptional configuration problem in 
the hairpin is even more serious than in the valence correlators. The MQA procedure allows 
the first detailed study of the time-dependence of the hairpin correlator. The most striking 
feature of this time-dependence is that it is quite accurately described at all time separations 
by a the simple chiral Lagrangian diagram consisting of two pion propagators on either side 
of a momentum independent mass insertion. Values obtained for the field renormalization 
parameter q;$, which parametrizes the momentum dependence of the hairpin insertion, 
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are very small and consistent with zero. Perhaps even more remarkable is the fact that the 
hairpin correlator exhibits an absence of excited state contamination, with the ground-state 
double-pole diagram giving a complete description of the correlator. This result is confirmed 
not only by the time-dependence of the hairpin correlator, but also by a detailed comparison 
of hairpin and valence correlators with both smeared and local sources, as discussed in Ref. 
[5. Since we know from the valence correlator that the local ip^^ip operator creates a state 
which includes a substantial excited state component along with the ground state pion, 
we conclude that these excited states are decoupled from the hairpin vertex itself. This 
may be viewed as a plausible extension of the OZI rule. Calculation of vector and axial- 
vector hairpin correlators |2J has confirmed that qq annihilation in these channels is highly 
suppressed compared to the anomaly-enhanced pseudoscalar hairpin diagram, as expected 
from OZI phenomenology. The absence of evidence for excited states in the pseudoscalar 
hairpin indicates that annihilation from these excited pseudoscalar states is similarly OZI 
suppressed. Only if the qq are in a Goldstone state do they have an unsuppressed pair- 
annihilation amplitude. The possibilities for extending these calculations to investigate 
the details of the connection between quark pair creation and annihilation and topological 
charge fluctuations are intriguing. 

A global fit to the pseudoscalar masses and decay constants obtained from the valence 
(nonsinglet) correlators as a function of the quark masses yields an estimate of the quenched 
chiral log parameter 6, the pion decay constant the slope parameter tq, and the Leutwyler 
parameters L5 and Lg (Table Oj). We have found that a chiral Lagrangian analysis at the 
one-loop level provides a good description of the data over the range of masses studied, for 
both equal and unequal quark mass. The fit values of the quenched chiral log parameter 
6, ~ .05 — .06, are rather small, but consistent with the values obtained directly from the 
hairpin diagram, averaged over the mass region studied. The chirally extrapolated value 
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of S (~ .10) evaluated from the hairpin correlator is somewhat larger, but still smaller 
than the phenomenological estimate of ~ 0.17. There may be some indication (cf. Fig. 
2) that the value of 5 is increasing for smaller lattice spacing so that, in the continuum 
limit, it might be closer to the phenomenological estimate. Similarly, the values of L5 (as) 
we obtain are somewhat larger than the recent results of Heitger et.al. |S]. In any case, 
since the lattice spacings we have studied are fairly coarse, the overall consistency of our 
fits suggests that, even for finite lattice spacing, the low energy dynamics is well-described 
by a chiral Lagrangian, with the main lattice spacing effects consisting of corrections to the 
parameters. 

The behavior of the nonsinglet scalar propagator at light quark mass provides a partic- 
ularly useful probe of chiral dynamics, and the success of the chiral Lagrangian description 
(Section 6) is impressive. The prominent r/'-vr loop contribution is completely determined 
in terms of the chiral parameters m^r, tq, and rriQ, extracted from the pseudoscalar analysis. 
The loop calculation agrees well with the lattice data in magnitude, time-dependence, and 
pion mass dependence. For the box size 2 fm) and pion masses we have used, the loop 
diagram exhibits a strong finite volume effect, and the agreement with the data is only 
satisfactory if the loop calculation is carried out in a finite volume. It would be interesting 
to carry out the numerical scalar propagator calculation on different size lattices to explic- 
itly observe this finite volume effect. In view of the agreement we see on a single size box 
over a wide range of pion masses, we would expect the chiral loop diagram to also give a 
good description of the finite volume dependence, as long as the box is still large compared 
to the QCD scale. The sensitivity of the r/'-vr loop diagram to finite volume effects makes 
it especially useful for probing the role of zero modes in finite volume calculations. For 
extremely large boxes (large compared to the chiral scale), zero modes and global topology 
should be irrelevant, e.g. ensemble averages over any fixed topology should converge to the 
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same result in the infinite volume limit. On the other hand, for volumes which are large 
with respect to the QCD scale but still comparable to the chiral scale, the exact zero modes 
contribute in an essential way to the dynamics described by the chiral Lagrangian. The 
techniques developed here can easily be used to study the role of zero modes and global 
topology in finite volume chiral dynamics. This can be done, for example, by studying the 
contributions to the ensemble average for various correlators as a function of the global 
topological charge. We have seen that the integrated anomaly method of Smit and Vink 
jl7j . after MQA improvement, provides a convenient way of estimating the global topologi- 
cal index v oi a, configuration. Calculation of u using exactly chiral (e.g. overlap or domain 
wall) fermions would be ideal, but the approach used here based on clover improved Wilson 
fermions is more economical and appears to be quite effective for studying issues which do 
not depend crucially on having exact integer valued z^'s. For example, the results for the 
topological susceptibility (Section 4) are in good agreement with other estimates. It would 
be interesting to carry out a more detailed investigation of the scalar and pseudoscalar 
correlators studied here, decomposing the ensemble averages into various topological charge 
sectors. This would provide useful insight into the role of global topology in chiral dynamics 
at finite volume. 

After properly accounting for quenched chiral effects in the scalar propagator, the oq 
masses (at Kc) are extracted. The results, rUaQ = 1.33(5) and 1.34(6) GeV for (3 = 5.9 and 
5.7 respectively are considerably larger than the observed ao(980) resonance mass. The 
value for lUaQ = 1.33(5) GeV suggests that there is a large effect when internal quark loops 
are included or that the observed ao(980) resonance is a distinct state, possibily a KK 
"molecule" (which would not appear in the quenched approximation), and not an ordinary 
qq meson. 
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